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Abstract. With the aid of stellar atmosphere models, we derive iron 
abundances [Fe/H] from the OGLE B, V, I photometry on the Mag- 
ellanic Cloud Cepheids. We show that in both clouds the average 
, metallicities of the first overtone variables are lower than those of the 

fundamental ones (by ~ 0.2 and ~ 0.3 dex in the LMC and SMC, 
respectively). Consequently, there is a correlation between the overall 
< I . [Fs/H] and luminosity; the lower luminosity stars tend to be also of 

^ ' lower metallicity. These metallicity dependencies are concordant with 

^ the ones derived for the two types of double-mode Cepheids from pul- 

sation theory. Yet another support of this dichotomy comes from the 
evolution theories that require lower metallicities for blue-looping low- 
luminosity stars than for high-luminosity ones. We also comment on 
the possibility of using period-luminosity-color relations to derive more 
accurate metallicities. 
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1 Introduction 



O . 

OO , It is well known that chemical composition plays a very important role in stellar 
' evolution and in various methods of stellar parameter computations, most remark- 
^ , ably in asteroseismological investigations. From the point of view of chemical com- 
' position, stellar clusters and certain populations of galaxies have been treated quite 
rS , often as chemically homogeneous ensembles. However, accumulating abundance 
' measurements show that there are chemical inhomogeneities even in globular clus- 
ters (e.g., Yong & Grundahl I2008P . where the uniform composition thought to be 
well-justified on the basis of common origin of the constituting objects. It has also 
been evident during the recent years ~ mostly due to the systematic spectroscopic 
surveys of Andrievsky and coworkers - that there is a considerable metallicity 
spread among the Galactic Cepheids and that the observed metallicity correlates 
with the distance from the Galactic center (Luck, Kovtyukh, & Andrievsky I2006( ) . 
The importance of the heavy element abundance in stellar pulsation has become 
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especially clear after 1992, with the verification of Simon's hypothesis on the iron 
opacity bump f Simon [T9821 Rogers & Iglesias [l992p . Perhaps the largest impact 
made by the opacity increase was the near elimination of the so-called 'Cepheid 
mass discrepancy' (Moskalik, Buchler & Marom ri992p and the excitation of the 
P Cephei stars (Moskalik & Dziembowski ll99^ . In addition, based on these revised 
opacities, there are theoretical results which show that the consideration of metal- 
licity spread is a must in various groups of Cepheids. Currently it has been shown 
that the observed periods of the first/second overtone double-mode Cepheids in 
the Magellanic Clouds can be fitted only if we assume that the metallicities of these 
low-luminosity stars are lower than their fundamental/first overtone counterparts 
(Kovacs l2006p . This observation supports the earlier result of Cordier, Goupil & 
Lebreton pOOSp . suggesting that Cepheids cannot 'blue-loop' at low- luminosities 
in the SMC, unless they have systematically lower metallicities than the higher 
luminosity ones. Stimulated by the above results, the goal of the present study is 
to use the photometric metallicities in the investigation of the metallicity spread 
among the Magellanic Cloud Cepheids. 

2 Metallicities from BVI colors 

Our method has already been described in Dekany et al. (|2008p and tested on 
a sample of 21 fundamental mode Galactic field RR Lyrae stars. Here we sum- 
marize the main steps of the method and present tests concerning the limits of 
applicability of the photometrically derived abundances. 

Assuming that the poorly known quantities, such as the turbulent velocity 
and convective mixing length are fixed to some generally acceptable value and the 
relative abundances are also given (e.g., scaled to the Sun), then stellar atmosphere 
models and basic pulsational relations yield the following set of equations: 

S-y = i^(Teff,5,[Fe/H]) , F -/ = G(Tcff, 5, [Fe/H]) , g = H{P) , (2.1) 

where the last equation is applicable in a wide range of parameters in the region 
of Cepheids (see Kovacs 2000). The functions F and G denote the dependence 
of the theoretical color indices on the physical parameters as given by the models 
of Castelli, Gratton & Kurucz (|1997p . All models have solar- type heavy element 
distribution, microturbulent velocity of 2 kms^^ and no convective overshooting. 
Quadratic interpolation among these models proved to be accurate enough to find 
the solution of set of Eqs. (2.1), which is searched for by minimizing the following 
quantity: 

= [A log Teff] ' + K-v A(B -!/)]'+ [w-i A(y -/)]', (2.2) 

where A log Tog is the difference between the log Tcff values computed from the 
B — V and V — I colors, whereas A{B — V) and A{V — I) denote the differences 
between the observed values and the ones given by the stellar atmosphere model. 
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The weights Wb-y and wy-i are set equal to 0.33 and 0.25, respectively, following 
the proportionality of log T^g in the approximate linear expressions of Kovacs & 
Walker p999p . The last two terms in V increases the stability of the method 
against observational noise. In general, we accept solutions only with T> < 0.001. 




[Fe/H]..t. [Fe/H], 



Fig. 1. Left: theoretical vs. computed photometric metallicities by shifting the theoret- 
ical V magnitudes by AV^ = —0.01 mag and use this (with the other input parameters 
left at their theoretical values) in the computation of [Fe/HJavi- Right: the same plot 
for the B color with AB = +0.01 mag. 

For variable stars it is a question of what kind of average magnitudes should 
we use if we need a good approximation of the static colors. We use the intensity 
averaged magnitudes mostly because of the easy availability of these values from 
the OGLE database. Apparently the intensity averages yield closer agreement 
with the static colors, at least for RR Lyrae stars, concerning B, V and K colors 
(Bono, Caputo & SteUingwerf IT995|) . 

The method was tested on various datasets. For a sample of 42 Pop. I dwarfs 
with spectroscopically determined [Fe/H], log 5 and Teff values and BVI photom- 
etry (e.g., Clementini et al. 119951 Tavlor [2003D . we got a = 0.14 dex between 
the spectroscopic and our photometric [Fe/H]. A much more limited sample of 10 
Galactic Cepheids yielded a = 0.19 dex. On the set of the 21 Galactic RR Lyrae 
stars mentioned above, we derived a = 0.15 dex. No apparent systematic dif- 
ferences were observed between the spectroscopic and photometric abundances. 
We conclude that the photometric abundances derived from the BVI colors are 
reasonably accurate on various types of stars. We will see that this assessment is 
confirmed also by the average metallicities derived for the MC Cepheids. 

To obtain additional information on the sensitivity of the photometric abun- 
dances, here we examine the differences obtained between the tabulated theoretical 
values and the ones obtained by various perturbations of the corresponding input 
parameters (as given by the stellar atmosphere models). 

In Fig. 1 we show two examples on the effect of changing the photometric zero 
points. We see a dependence on the metallicity both in the overall difference and 
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Table 1. Systematic and random effects for the photometric [Fe/H] and TaS 
Systematic errors: 



Parameter 


shift 


A [Fe/H] 


ATeff 


AV 


+0.01 


-0.31 + 0.18 


-34+12 




-0.01 


+0.41 + 0.28 


+31 + 12 


AB 


+0.01 


+0.25 + 0.21 


- 3 + 


7 




-0.01 


-0.17+0.12 


+ + 


6 


AI 


+0.01 


+0.22 + 0.15 


+33 + 


5 




-0.01 


-0.18 + 0.12 


-34 + 


6 


A log 5 


+0.10 


-0.01 + 0.03 


+ 8 + 


7 




-0.10 


+0.01 + 0.04 


- 8 + 


8 


Random errors: 










Parameter 


a 


CT[Fe/H] 


o-(T'off) 


a{B -V)^ ct{V - I) 


0.01 


0.26 


35 




a[E{B - V)] 


0.01 


0.05 


44 





Notes: The table shows the errors (computed minus model values) of the photometric [Fe/H] 
and Tefj for various assumed errors in the input parameters. All other input parameters are 
left at their theoretical values, except for the parameter indicated. The symbol A in the 
parameter column means 'shifted minus theoretical' values. We used Gaussian noise in the 
random simulations (last two rows). In computing the assemble values, we used 4cr clipping. 
The errors depend on the actual parameter regime (see, e.g.. Fig. 1). The errors shown in this 
table are average values including the full parameter range. 

also in the size of scatter (due to the various values of Tcfi and log g). On the other 
hand, it is conceivable that if we are interested in relative abundances, then the 
method is more accurate, especially at higher metallicities. 

To get a closeup on the effects of other possible systematic shifts and random 
errors, we summarize the resulting [Fe/H] errors in Table 1 (for completeness, 
we also show Toff). As it is expected, photometric errors play the most important 
role, while errors due to reddening are smaller. Gravity estimates can be especially 
erroneous without affecting the derived metallicity. 

3 [Fe/H] distributions in the MCs 

We applied the above method of [Fe/H] determination to the Cepheid database 
published by Udalski et al. (|1999p . We used their periods, intensity-averaged 
magnitudes and reddening values as given in their *.tah filesQ d To avoid large 
errors due to observational noise, as mentioned above, we accepted solutions only 
if they satisfied the T) < 0.001 criterion. In addition, we also employed a 3a 
clipping to filter out extreme values. The resulting [Fe/H] distributions are shown 
in Fig. 2. 



^ ftp: // sirius.astrouw.edu. pi / ogle / ogle2 / var_stars /Imc / cep / catalog 
^ftp:/ /sirius. astrouw.edu.pl/ogle/ogle2/var_stars/smc/cep/catalog/ 
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[Fe/H] 



[Fe/H] 



Fig. 2. Metallicity distributions computed directly from the observed BVI colors. 



Table 2. Average photometric metallicities for the Magellanic Cloud Cepheids 



Cloud 




[Fe/H]Fu 


Nfo 


[Fe/H]Fo 


[Fe/H]Fu-[Fe/H]Fo 


LMC 


676 


-0.27 ±0.02 


437 


-0.43 ±0.02 


0.16 ±0.03 


SMC 


941 


-0.97 ±0.02 


519 


-1.23 ±0.03 


0.26 ±0.03 



It is clearly seen that there is a surplus of first overtone mode variables at the 
low- metallicity end of the distributions in both clouds. The same is true for the 
fundamental mode variables at the high-metallicity end. In both clouds the metal- 
licity spread is large but it is clear that most of this comes from observational 
errors (e.g., the spread is larger in the SMC, because of the larger photometric 
errors). When computing the average metallicities, we get the result shown in Ta- 
ble 2. Based on the formal errors of the mean values, the differences in the overall 
metallicities seem to be significant and are exhibited in the same sense in both 
clouds. We note that: (i) the average photometric metallicities are quite close to 
the 'canonical' values (e.g., Udalski et al. [1999D; (ii) the derived lower metallicity 
for the SMC FO variables is the same as suggested by Cordier et al. (|2003p for 
the low-luminosity Cepheids. 



4 Combination with the PLC relations? 

In the direct computation of [Fe/H] from the observed color indices, large errors are 
introduced due to the high sensitivity to the observational errors. In principle, a 
considerable fraction of this error would be possible to eliminate if we could employ 
the period-luminosity-color (PLC) relations and use only the /-band magnitudes 
and the period to compute the color indices. Let us assume that for each reddening- 
free magnitude we have the same type of relations for all Cepheids pulsating, e.g., 
in the fundamental mode: 



Wb-y 



= /-1.86(B-y) = ao + ailogPo , 



(4.1) 
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Wy-i = J-1.50(F-/) = 6o + 6ilogPo , (4.2) 

where the and hi coefficients are independent of the metallicity. Similar relations 
are assumed to be hold also for the FO variables. The obvious advantage of the 
above relations is that the errors in the computed color indices become highly 
correlated through the errors of the directly used color /. Indeed, if we test this 
method on the stellar atmosphere models of Castelli et al. (|1997p . then, with 
an error of ct/ = 0.03, we get cr[p(,/H] = 0.06 (we recall that independent errors 
of 0.01 in the color indices cause errors of cr [fo/h] ~ 0.26 — see Table 1). We 
note in passing that Tgff will not become better estimated. With the above noise 
we get tJT^ft = 65. This high insensitivity of [Fe/H] against observational noise 
can be understood by the inspection of the linear approximation of the color- 
temperature dependence in the above colors. The coefficients of the color terms 
are nearly cancelled out when solved for the metallicity with the aid of the above 
PLC relations (see Kovacs & Walker [19991) ■ 

Although the above method is very promising, there are considerable difficulties 
when looking at the assumptions in more detail. First of all, as it follows from the 
M, L and Tofi dependence of the pulsation period, the PLC relations should have 
some internal scatter, because the mass-luminosity relation is not a strict one. 
Therefore, the derived relation may depend on the stellar population studied and 
cannot be employed in a straightforward way to other systems. The population 
dependence is obvious in the case of PL relations and is widely discussed in the 
literature. Although PLC relations have been shown to be much less sensitive to 
such effects, in fine details we might face with similar problems as in the case of the 
PL relations. This question has recently been studied by Fiorentino et al. (|2007p . 
They showed that the PLCs based on the B, V colors have significant systematic 
dependence on the metallicity. The effect is smaller but not negligible also in 
the near infrared colors. This has led the authors to estimate the metallicities of 
Galactic Cepheids which they showed to exhibit a fine correlation with the direct 
spectroscopic data of Andrievsky and coworkers. 

There are also ambiguities in the derivation of the PLC relations. Udalski et 
al. (|1999p employ a 2.5(T-clipping and introduce various cuts in their MC sample. 
In this way they end up with the same slope in WV-i for the FU variables in both 
clouds. However, they note that the slopes for the FO variables might be different. 
Also, they did not study Wb-v, presumably due to the paucity of the B data. 

We derived the PLC relations in both unrcddened colors as given by Eqs. (4.1) 
and (4.2). We do not apply any period constraint on the sample. For each class of 
variables the first 40 strongest outliers have been omitted iteratively. Table 3 gives 
the resulting regression parameters. By comparing the same types of variables 
in the two clouds, we see that there is a significant difference in hi for the FU 
variables. Although we use a slightly different expression for Wy-i than Udalski 
et al. (|1999p . we think that this result is mainly due to the drop of the period 
constraint. We recall that Udalski et al. (|1999p derived the same slope in both 
clouds within fairly narrow error limits. On the other hand, we agree with them 
on the FO variables: the SMC variables have a steeper slope in Wy-i than the 
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Table 3. PLC relations for the Magellanic Clouds Cepheids 







I -1M{B -V) 


= ao + fli logP 






Cloud 


Mode 


ao 


ai 




N 


LMC 


FU 


16.021 ±0.013 


-3.590 ±0.019 


0.109 


732 


SMC 


FU 


16.677 ±0.008 


-3.590 ±0.016 


0.166 


1211 


LMC 


FO 


15.444 ±0.008 


-3.495 ±0.023 


0.094 


469 


SMC 


FO 


16.150 ±0.007 


-3.604 ±0.028 


0.150 


734 






I -1.50{V-I) 


= 6o±6i logP 






Cloud 


Mode 


bo 


bi 




N 


LMC 


FU 


15.915 ±0.007 


-3.314±0.011 


0.064 


732 


SMC 


FU 


16.521 ±0.006 


-3.419 ±0.012 


0.132 


1211 


LMC 


FO 


15.405 ±0.006 


-3.397 ±0.015 


0.063 


469 


SMC 


FO 


16.011 ±0.006 


-3.556 ±0.025 


0.134 


734 



ones in the LMC. A similar effect is seen for Wb-y- Accidentally, the FU variables 
behave in the same way in both clouds in respect of the above parameter. It is 
also important to note that the Wb-v index yields ^ 0.15 mag larger value for 
the relative distance modulus of the clouds than the 'canonical' value of 0.51 mag 
derived from the WV-i index. 

If we disregard all the above inconsistencies and use the PLC method to com- 
pute [Fe/H], we get results confirming the metallicity difference between the FO 
and FU Cepheids in both clouds. For the LMC we get average metallicities of 
—0.22 and —0.36 for the FU and FO variables, respectively. For the SMC we 
get —1.24 and —1.63, that are considerable lower than the ones obtained directly 
from the color indices (however, the relative difference between the FU and FO 
variables remains nearly the same). In both clouds (but especially in SMC) we 
get a considerable decrease in the total metallicity range. 

5 Conclusions 

To lend further support to the hypothesis of systematic metallicity spread among 
the Magellanic Cloud Cepheids, we attempted to derive photometric metallicities 
from the OGLE BVI photometry. Aided by the stellar atmosphere models of 
Castelli et al. (|1997p we have shown that first overtone variables have systemati- 
cally lower abundances than fundamental ones. The effect is larger in the SMC, 
amounting to ~ 0.3 dex (vs. ~ 0.2 dex in the LMC). This result is in agreement 
with the ones suggested by the evolutionary models of Cordier et al. (2003) on the 
condition of instability strip crossing of low- luminosity Cepheids in the SMC. Also, 
the discrepancy between the pulsational and evolutionary ML relations (Beaulieu 
et al. I200ip can be largely cured by assuming lower metallicities. Yet another, 
independent result in fitting the periods of the SO /FO Cepheids in the Magellanic 
Clouds shows that they require lower metallicities than their FO /FU counterparts 
(Kovacs G. I2006p . We think that all these results point toward the reality of the 
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metallicity dichotomy between the FO and FU variables. Obviously, a massive 
spectroscopic survey of the Magellanic Clouds Cepheids (similar to that of Mot- 
tini et al. 2006.) would be of great value, because the photometric abundances 
are of low accuracy and cannot tell the degree of the dichotomy and its more 
precise relation to the luminosity. We also attempted to use the PLC relations 
to derive more accurate color indices and thereby more precise abundances. Un- 
fortunately this method remains somewhat doubtful due to controversies in the 
observed relations and due to inherent scatter of the theoretical ones (Fiorentino 
et al. [2007|) . 
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